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Abstract. The equations of the relativistic causal Newton gravity law for the planets of the 
solar system are studied in the approximation when the Sun rests at the coordinates origin 
and the planets do not interact between each other. 



1 Introduction 



The Newton gravity law requires the instant propagation of the force action. The special 
relativity requires that the propagation speed does not exceed the speed of light. If the 
propagation speed is independent of the gravitating body speed, then it is equal to that of 
light. The special relativity requires also the gravity law covariance under Lorentz transfor- 
mations. Poincare [1] tried to find such a modification of the Newton gravity law. (Poincare 
considered two mathematical problems in XX century as principal: "to create the math- 
CSj ; ematical basis for the quantum physics and for the relativity theory") The gravity forces 

of two physical points should depend not on its simultaneous positions and speeds but on 
£T) ■ the positions and the speeds at the time moments which differ from each other in the time 

interval needed for light covering the distance between the physical points. The gravity force 
acting on one physical point may depend also on the acceleration of another physical point 
at the delayed time moment. The relativistic Newton gravity law was proposed in the paper 
[2]. This law for the two physical points has the form 



d ((, .-2 <&* 2 V 1/2 ix k\ _«» ^ _-i<&* 




dt 



(U , '/"" E c l ^F j; , u (x k , Xj ), (1.1) 



j, k — 1, 2, j k, fi = 0, 3. The world line x^it) satisfies the condition x° k {t) = ct; c is the 
speed of light; the diagonal 4x4- matrix i] fiu = i]^ u , rf° = —rj 11 = —rf 2 = —t] 33 = 1; the 
strength Fj.^^XkyXj) is expressed through the vector potential 

j-, / \ 9A.j- v (^Xk, Xj) dAj-^^XkjXj) , . 

fj;fiv{Xk,Xj) — ~3x^ dx~l ' 

Aj- ttx {x k , Xj) = rj^mjG (c|x* - x,(t')| - - x )(t'))^x){t')^ , (1.3) 



*This work was supported in part by the Russian Foundation for Basic Research (Grant No. 12-01- 
00094), the Program for Supporting Leading Scientific Schools (Grant No. 4612.2012.1) and the RAS 
Program "Fundamental Problems of Nonlinear Mechanics." 
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* , = c" 1 (a;2-|x fc -x J -(0|), j, k — 1,2, j^k; 

the gravitation constant G = (6.673 ± 0.003) • 10 _11 m 3 A;g _1 s~ 2 and rrij is the j body mass. 
For a resting body world line (xj(t) = ct and the vector Xj(t) is constant) the vector potential 
(1.3) coincides with the Coulomb vector potential 

Aj. (x k ,Xj) = mjGl^k - Xj(c -1 a;2)r\ Aj.i(x k ,Xj) = 0, i = 1,2,3. (1.4) 

If the velocities of bodies are small enough to neglect their squares compared with the square 
of the light speed and it is possible to neglect also the time interval c _1 |x fc — x.j(t')\, then 
the vector potential (1.3) is nearly equal to the Coulomb vector potential (1.4). The vector 
potential (1.3) was proposed by Lienard (1898) and Wiechert (1900) as the generalization of 
the Coulomb vector potential (1.4). The substitution of the Coulomb vector potential (1.4) 
into the right-hand side of the equation (1.1) for fi = 1,2,3 yields the right-hand side of 
the Newton gravity law equations. The equation (1.1) multiplied by (1 — c~ 2 \dx.k/dt\ 2 ) 
transforms as the vector. The equations (1.1), (1.2) with the Lienard - Wiechert vector 
potential (1.3) are the relativistic version of the Newton gravity law equations. 

Sommerfeld ([3], Sec. 38): "The question may arise: what is the relativistic form of the 
Newton gravity law? If the law is supposed to have a vector form, this question is wrong. 
The gravitational field is not a vector field. It has the incomparably complicated tensor 
structure." The Newton gravity law equations and the equations (1.1) - (1.3) define the 
interactions. The body interacts only with another body. If two bodies create the common 
gravitational field with the vector potential A^{x, Xi)+A 2; ^(x, x 2 ), any body should interact 
with itself and we obtain the infinity in the equations (1.3), (1.4) at x k = Xj. The notion of 
gravitational field with the vector potential A\.^(x, x\) + A 2 - lJL {x, x 2 ) is not compatible with 
the Newton gravity law and with the relativistic Newton gravity law (1.1) - (1.3). 

The delay c _1 |x fc — Xj(t')| in the relation (1.3) provides the causality condition according 
to which some event in the system can influence the evolution of the system in the future only 
and can not influence the behavior of the system in the past, in the time preceding the given 
event. The delay c _1 |xfc — Xj(£')| in the relation (1.3) is very important: one celestial body 
is a good distance off another celestial body. Poincare [1]: "It turned out to be necessary to 
consider this hypothesis more attentively and to study the changes it makes in the gravity 
laws in particular. First, it obviously enables us to suppose that the gravity forces propagate 
not instantly but at the speed of light." The general relativity does not take into account 
the causality condition and the delay. 

99.87% of the total mass of the solar system belongs to the Sun. We consider the 
relativistic causal Newton gravity law equations [2] for the planets of the solar system in the 
natural approximation when the Sun rests at the coordinates origin and the planets do not 
interact between each other. In this approximation the problem of planet relativistic motion 
was solved in the paper [2]. The planet orbits were given by the formulas which differ 
from the formulas defining the ellipses in the precession coefficients only. The precession 
coefficients for the solar system planets are practically equal to one. The similar orbits with 
another precession coefficients are considered in the general relativity ([4], Chap. 40, Sec. 
40.5, relations (40.17), (40.18)). In the beginning of the XVII century Johannes Kepler by 
making use of Tycho Brahe (1546 - 1601) astronomical observations found that the planet 
orbits are elliptic in the coordinate system where the Sun rests (Nicolaus Copernicus (1543)). 
The intensive astronomic observations from the middle of the XIX century and the radio- 
location after 1966 discovered the advances of orbit perihelion for different planets. 
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In the general relativity the observed value for the Mercury's perihelion advance is ob- 
tained by means of addition the advance of Mercury's perihelion ([4], Chap. 40, Sec. 40.5, 
Appendix 40.3) calculated in the Newton gravity theory and the advance of Mercury's per- 
ihelion calculated for the orbits ([4], Chap. 40, Sec. 40.5, relations (40.17), (40.18)). The 
orbits ([4], Chap. 40, Sec. 40.5, relations (40.17), (40.18)) are the approximate solutions of 
the geodesic equation for the chosen metrics ([4], Chap. 40, Sec. 40.1, relation (40.3)). It is 
not obvious that we can add the advance of Mercury's perihelion obtained for the orbits ([4], 
Chap. 40, Sec. 40.5, Appendix 40.3) and for the orbits ([4], Chap. 40, Sec. 40.5, relations 
(40.17), (40.18)). It seems natural to obtain the advance of Mercury's perihelion, observed 
from the Earth, by making use of the Mercury and Earth orbits ([4], Chap. 40, Sec. 40.5, 
relations (40.17), (40.18)) calculated without Newton gravity theory. In order to calculate 
the advance of Mercury's perihelion we need to know also the time dependence of the orbit 
([4], Chap. 40, Sec. 40.5, relations (40.17), (40.18)) radius. In this paper we study the 
similar orbits of the relativistic causal Newton gravity law [2]. We shall show in this paper 
that the value of the Mercury's perihelion advance, observed from the Earth, depends on 
the perihelion angles of the Mercury and Earth orbits. The perihelion angle of the planet's 
orbit depends on the planet's perihelion point due to the precession coefficient in the planet's 
orbit formula. For the experimental verification of the equations (1.1) - (1-3) the perihelion 
angles of the Mercury and Earth orbits are needed. 



2 Causal Coulomb and Newton laws 

The relativistic Lagrange law is the particular case of the relativistic Newton second law 

dt d ( dt dx»\ ^ dtdx ai dt dx ak n 

mc-- - — — + gcT 1 V V ?? w F„ Ql ... Q . (x)- — = 0, 2.1 

dsdt\ds dt y t n„ „ n d s dt ds dt ' v ; 

dt / 9 ,9\- 1 /2 ,■ dx l . „ n „ 
_ = (cMv p) ,„ = _,, = 1,2,3. 

where /i — 0, 3 and the world line x^(t) satisfies the condition: x°(t) = ct. The force is the 
polynomial of the speed in the equation (2.1). It is necessary to define the series convergence 
for the force as an infinite series of the speed. The second relation (2.1) implies the identities 

Jt, f dt dx a \ 2 ^ dtdx a dt d ( dt dx a \ 

ho ^ \d~s~dTj ~ ' h^dslttdsJt [d~s^r) { } 



The equation (2.1) and the second identity (2.2) imply 

dtdx ai 

w - afe + l(a %^T' ds dt 



* . ,dtdx ai dtdx a ^ 

E E ^W^-- — = «■ ( 2 - 3 ) 

fe = ai,...,a k+1 =0 



Let the functions F ai ... ak+1 (x) satisfy the equation (2.3). Then three equations (2.1) for 
/i — 1, 2, 3 are independent 

j n 
m ± ((i _ C - 2 | V | 2 )-V2^ _ qc -i (c 2 - M 2 )-^- 1 )/ 2 
dt k = o 

( A , ,dx ai o 

x 



E F iai ... ak (x)^-.-.^-\=0, i = 1,2,3. (2.4) 



dt dt 

y ai,...,a k = U 
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The following lemma is proved in the paper [2]. 

Lemma. Let there exist a Lagrange function L(x, v, t) snc/i t/iat for any world line x^{t), 
x°(t) = ct, the relation 



d dL 


dL 


didv* ' 


dx i 






(c 2 - 


|v| 2 )- 



m— 

dt 



A N 
" ''I - c-VlVV) - qc' 1 £ 

fc = 
ai,...,a)t =U 



/10/ds /or any i = 1, 2, 3. Then the Lagrange function has the form 
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L(x, v, i) = -mc 2 (l - c~ 2 |v| 2 ) 1/2 + Q M x , *)c _ V + gA(x, t) (2.6) 

i = i 

and the coefficients F iai ... ak (x) in the equations (2.4) are 

F iai ... ak (x) = 0, k ^ 1, % = 1,2,3, «i, ...,a fc = 0, ...,3, (2.7) 

&4 3 -(x,t) ^(x,t) . 
^W" dx i dxi , M-M,*, 

dA (x,t) l dJU^t) 
F i0 (x) - - i-l,2,3. (2.8) 

We define the coefficients 

F 00 (x) = 0, F 0l (x) = -F i0 (x), i = l,2, 3. (2.9) 

Then the identity 

A dtdx a dt dx 13 

? F ^ x) dSlTTslir = (2 - 10) 

a,p = 

of the type (2.3) holds. By making use of the second identity (2.2) and the relations (2.8) - 
(2.10) we can rewrite the equation (2.4) with the coefficients (2.7), (2.8) as the relativistic 
Newton second law with Lorentz force 

dt d f dt dx» 

mc 



ds dt 



(dtdx»\ „„ 3 ,dtdx u 



^<») = ^-^U, = 0.....S. (2,1) 

For the relativistic Lagrange law the interaction is defined by the product of the charge q 
and the external vector potential A M (x, t). 

Let a distribution e (x) e S"(R 4 ) with support in the closed upper light cone be a 
fundamental solution of the wave equation 

-(d x ,d x )e (x) = S(x), (d x ,d x ) = (J^) - £ (^) • (2-12) 

We prove the uniqueness of the equation (2.12) solution in the class of distributions with 
supports in the closed upper light cone. Let the equation (2.12) have two solutions e^\x), 
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e'o'\x). Since its supports lie in the closed upper light cone, the convolution is defined. Now 
the convolution commutativity 

J d 4 xd*ye£\x - y)ei 1 \y) ( j ) (x) = J d A xd A ye^ \x)e^\y)<t>{x + y) (2.13) 

implies these distributions coincidence: 

e { i\x) = -(d x ,d x ) J d^\x-y)e i i\y) ) (2.14) 

j, k — 1, 2, j 7^ k. Due to the book ([5], Sect. 30) this unique causal distribution is 

e (x) = -(2n)- 1 6(x°)6((x,x)), (2.15) 



(s,y) = *y-E*V, W = ln 
k = i l u > 



x > 0, 
x < 0. 



The relativistic causal Coulomb law is given by the equations of the type (2.11) 

m k j t ( (l -c- 2 ^ 2 Y' 2d 4\ = -*«T E c-^F,^(x fc ,x,), (2.16) 



j, k = 1,2, j 7^ k, where the strength Fj-^(xk,Xj) is given by the relation (1.2) with the 
Lienard - Wiechert vector potential of the type (1.3) 



3 

Aj. tli (x k ,Xj) = -inqjK ^ r)^ J dte (x k - Xj(t)) 



dx]{t) 
dt 



u = 

1 



- qj K V „ (c|x fc -x,(t)i - E(4 



(2.17) 



t=t(o) 



4-ct(0) = |x fe -x,(t(0))|. 

Here X is the constant of the causal electromagnetic interaction for two particles with the 
charges qj. The support of the distribution (2.15) lies in the upper light cone boundary. The 
interaction speed is equal to that of light. It is easy to prove the second relation (2.17) by 
making change of the integration variable 

x°-ct(r) = (|x fc -x,(t(r))| 2 + r) 1 /2. ( 2 .1 8 ) 

For r = the relation (2.18) coincides with the third relation (2.17). 

The equations (2.16), (1.2), (2.17) are the relativistic causal version of the Coulomb law. 
The Lorentz invariant distribution (2.15) defines the delay. The Lorentz invariant solutions 
of the equation (2.12) are described in the paper [2]. By making use of these solutions 
it is possible to describe the Lorentz covariant equations of the type (2.16), (1.2), (2.17). 
The equations (2.16), (1.2), (2.17) are Lorentz covariant and causal due to the distribution 
(2.15). The quantum version of the equations (2.16), (1.2), (2.17) is defined in the paper [6]. 
The solutions of these causal equations do not contain the diverging integrals similar to the 
diverging integrals of the quantum electrodynamics. 
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For a world line Xj(t) we define the vector proportional to — r]^(d x , d x )Aj. fl {x, Xj) 

r dx^it) r dx^(t) 

J»(x, Xj ) = -(d x ,d x ) j dteo(x - Xj (t))-^ = J dt5(x- Xj (t))—±l = 

( C_V )) 6 ( X " Xj ( C ~ 1X °)) ' = °' 3 - (2 ' 19) 
The condition x^(t) = ct implies the continuity equation 

9 fA*. x,) = - £ (±x) (c- V)) ±6 (x - x, (c-4)) , 



dx° 

_d_ 

dx 



jj^xj) = (jLx) (c-^fj ^(x-x,- (c-V)) , i = 1,2,3, (2.20) 

x . ) = . (2.21) 



The integration of the relation 

e (x - xtf))-^ = J d*ye (x - y)S(y - ^(t))^ 1 (2-22) 
along the world line Xj (t) yields 

dx^(t) r 

dte (x - ^-(t))-^ = J d 4 ye (x - y)r(y, Xj ). (2.23) 



The relations (2.21), (2.23) imply the gauge condition for the vector potential (2.17) 

dx^ 



Due to the gauge condition (2.24) the tensor (1.2), (2.17) satisfies Maxwell equations with 
the current proportional to the current (2.19). 

The substitution K = —G and two positive or two negative gravitational masses q\ = 
±toi, q2 = ±m2 into the equations (2.16), (1.2), (2.17) yields the relativistic causal Newton 
gravity law (1.1) - (1.3). By changing the constants K = —G, qi = ±mi, q 2 = ±m 2 
in the equations from the paper [6] we have the quantum version of the equations (1.1) - 
(1.3). The substitution K = —G and also one positive and one negative gravitational masses 
q\ = ±777i, q2 = T m 2 into the equations (2.16), (1.2), (2.17) yields the galaxies scattering 
with an acceleration. Einstein [7]: "The theoretical physicists studying the problems of the 
general relativity can hardly doubt now that the gravitational and electromagnetic fields 
should have the same nature." 



3 Advance of Mercury's perihelion 



Due to the paper [2] the relativistic causal Newton gravity law for the solar system has the 
form 



dx k 



dt 



^ dx? 



dt 



3 A-v 

-f Ec" 1 ^ E l-^ir, .,-,). (3.1) 



u = 



J I '•»•./--/■• 



We give the number k — 1 for Mercury, the number k = 2 for Venus, the number k = 3 for 
the Earth, the number k — 4 for Mars, the number k = 5 for Jupiter, the number k = 6 for 
Saturn, the number k = 7 for Uranus, the number k = 8 for Neptune, the number A; = 9 for 
Pluto and the number k — 10 for the Sun. 

99.87% of the total mass of the solar system belongs to the Sun. We consider the Sun 
resting at the coordinates origin (Nicolaus Copernicus (1543)). Substituting the Sun world 
line x%(t) = ct, x\ (t) = 0, i = 1,2,3, into the equalities (1.2), (1.3) we have 



FiO;ij(x; x 10 ) =0, i,j = 1,2,3, F 10 . i0 (x; x w ) = -m 10 G|x| V, i = 1,2,3. 



(3.2) 



Substituting the Sun world line x%(t) = ct, x\ (t) = 0, % = 1, 2, 3, and the elliptic orbits into 
the expressions (1.2) and (1.3) it is possible to show that the values of strengths F 10 . i0 (x k ; x 10 ) 
considerably exceed the values of strengths Fj. iu (x k ,Xj) for any k,j = 1,...,9, k ^ j. It is 
possible to show also that the values of strengths Fj^ v (x\Q\ Xj) are negligible. We neglect 
the action of any planet on all of the other planets and the Sun. Then the Sun rests at the 
coordinates origin. Due to the relations (3.2) in the coordinates system where the Sun rests 
at the coordinates origin the first nine equations (3.1) have the form 




1 -c 



dx k 



dt 



-1/2 



dx\ 
~dl 



-m 10 G\x k \ 



1,2,3, k = 1,...,9 



(3.3) 



It is shown in the paper [2] that the following values 



Mi(x fe ) = e m 1 •' k 
»J=1 



, i dx{ _ j dx\ 

ijl I Xi. - - ' i 



dt 



dt 



R 






dt 



-1/2 



B 


dx k 




dt 


- mioG 





-1/2 



, / = 1,2,3, 



(3.4) 



(3.5) 



are conserved for the equations (3.3). The antisymmetric in all indices tensor e^i has the 
normalization £123 = 1. The conservation of the vector (3.4) is the relativistic second Kepler 
law. The vector x& is orthogonal to the constant vector (3.4). We introduce the polar 
coordinates in the plane orthogonal to the vector (3.4) 



x]_. k (t) = r k (t) cos (j) k (t), x\. k {t) = r k (t) sm<f) k (t), k = l,...,9. 
Let the constants (3.4), (3.5) satisfy the inequalities 

c 2 \M(x k )\ 2 -mi G 2 >0, 
|M(x fc )| 2 ((£(x fc )) 2 - c 4 ) + mj G 2 c 2 > 0, k = 1, 9. 
Due to the paper [2] the equations (3.3) have the solutions 



Pk 
r k (t) 



1 + e fc cos(7 fc (0 fc (t) - <j> k . fl )) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



where (fi k; o is the constant perihelion angle and the constants 

Pk = (c 2 |M(x fe )| 2 - m 2 10 G 2 )(m 10 GE(x k )) 



-1 



e k = (c 2 |M(x fe )| 2 ((£(x fe )) 2 - c 4 ) + m 2 l0 G 2 cy/ 2 (m 10 GE(x k ))-\ 
lk = (c 2 |M(x fc )| 2 - m 2 G 2 ) 1 / 2 (c|M(x fe )|)- 1 , k = 1, 9. 



(3.10) 
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The equations (3.9) are the relativistic first Kepler law. The orbit (3.9) is not periodic in 
general. The substitution of the vector (3.6) with r k (t) = a k , <pk{t) = Uk{t — t k -o) in the 
equations (3.3) yields the relativistic third Kepler law: 

(1 - c- 2 alul)- l ' 2 alujl = m w G. (3.11) 

The equations (3.9) define the trajectory of motion but do not define the time dependence of 
this trajectory. Let the constants (3.4), (3.5) satisfy the inequality (3.8) and the inequalities 

(£(x fe )) 2 < c\ fc = l,..,9. (3.12) 
Due to the paper [2] the equations (3.3) have the solutions with the constant parameter 

r fc (60 = m 10 G£(x fe )(c 4 - (£(x fc )) 2 )^ 
x(l + e fc sin((c 4 -( J E(x fc )) 2 ) 1 /V 1 a)), 
= m lQ G{E^k)fc-\c A - (£(x fe )) 2 )- 3/2 
x (c 3 (£(x fc ))- 2 (c 4 - (£( Xfc )) 2 ) 1/2 fe - &;<,) 

-e fc cos((c 4 - (E(x k )) 2 ) 1/2 c-%)), k = 1,...,9. (3.13) 

Let us express the constants in the equations (3.9), (3.13) trough the astronomical orbit 
data. The orbit eccentricities: e 1 = 0.21, e 2 = 0.007, e 3 = 0.017, e 4 = 0.093, e 5 = 0.048, 
e 6 = 0.056, e 7 = 0.047, e 8 = 0.009, e 9 = 0.249. Therefore < e k < 1, k = 1,...,9. Let us 
suppose E(x k ) > 0, k = 1, ...,9. The curve (3.9) is an ellipse with a precession. The focus 
of this ellipse is the coordinates origin. The major and minor "semi - axes" are equal to 

a k = Pk (l - el)- 1 = m 10 GE(x k )(c 4 - (£(x fc )) 2 )-\ (3.14) 

b k = a k (l - el) 1 / 2 = (c 2 |M(x fe )| 2 - m^G 2 ) 1 / 2 ^ 4 - (^(x,)) 2 )- 1 / 2 , k = 1, 9. (3.15) 

The inequalities (3.7) and e 2 < 1 imply the inequality (3.12). Hence the Eqs. (3.13) hold. 
For the parameters £ fc;± = ±(7r/2)c(c 4 — (_E(x fc )) 2 ) -1 / 2 we have the extremal radii 

r fc (& ; ±) = m 10 GE(x fc )(c 4 - (E^)) 2 )- 1 ^ ± e k ), k = 1, 9. (3.16) 

Hence, the "period" of the motion along the ellipse (3.9) is equal to 

T k = 2\t k (^ k . + ) - i fc (&._)| = 2nm 10 Gc 3 (c 4 - (E(x fe )) 2 )- 3 / 2 , k = 1, 9. (3.17) 

Let us define the mean "angular frequency" Uk = lnT^ 1 . The relation (3.17) implies 

u k = (c 4 -(E^k)) 2 ) 3/2 (rn w Gc 3 )-\ 
(E(x k )) 2 = c 2 (c 2 - (w fc m 10 G) 2/3 ), fc = l,...,9. (3.18) 

The substitution of the expression (3.18) into the equality (3.14) yields 

m 10 G = u 2 k al(2- 1 (l + a k (l-(2a k u k c- 1 ) 2 ) 1 / 2 )y 3/2 , a k = ±1, k = l,...,9. (3.19) 

Let c — >■ oo. Then mioG = w|a|(2 _1 (l + Ofc))~ 3 / 2 . For a k — 1 this expression agrees with 
the third Kepler law m w G = oj 2 a\. Choosing a k — 1 in the relation (3.19) we get the 
"relativistic third Kepler law" for the orbit (3.9) 

u k 2 afm w G= (2- 1 (l + (l-4o; 2 a 2 c- 2 ) 1 / 2 )) _3/2 ^l + ^ 2 a 2 c- 2 , fc = l,..,9. (3.20) 
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According to the book ([4], Chap. 25, Sec. 25.1, Appendix 25.1) the values u) k a\c~ 2 = 1477m 
for k = 1,2,3,4,6 (Mercury, Venus, the Earth, Mars and Saturn), the values ufafc~ 2 = 
1478m for 1 — 5,8 (Jupiter and Neptune), the value u 2 ajC~ 2 = 1476m for Uranus, the value 
w|agC" 2 = 1469m for Pluto; the major semi-axes a x = 0.5791 ■ 10 n m, a 2 = 1.0821 • 10 n m, 
a 3 = 1.4960 • 10 n m, a 4 = 2.2794 • 10 n m, a 5 = 7.783 • 10 n m, a 6 = 14.27 • 10 n m, a 7 = 
28.69 • 10 n m, a 8 = 44.98 • 10 n m, a 9 = 59.00 • 10 n m. The values co 2 k a 2 k c~ 2 = a k l ■ ula 3 k c~ 2 , 
k — 1, 9, are negligible and therefore the Sun mass values (3.20) obtained in the relativistic 
Kepler problem good agrees with values 0J k a\ obtained in Kepler problem. 
The substitution of the expression (3.20) into the equality (3.18) yields 

c" 4 (£(x fc )) 2 = 1 - 2u 2 k a 2 k c- 2 (l + (1 - 4 W 2 a 2 c" 2 ) 1 / 2 )" 1 » 1 - c^c" 2 , k = 1, 9. (3.21) 
By making use of the relations (3.10), (3.14), (3.15), (3.20), (3.21) we have 

7, = (l + 4^a 2 c~ 2 (l - e 2 )- 1 (l + (1 - ^c" 2 ) 1 / 2 )") - 

1 - 2'W k a 2 k c- 2 (l - e 2 )" 1 , k — 1, 9. (3.22) 

The value 2~ 1 a; 2 .a 2 ,c~ 2 (l — e 2 ) -1 pa 1 — ^ k is maximal for Mercury: 1 — ^ pa 1.3341 ■ 10~ 8 . 
The precession coefficients (3.22) of the orbits (3.9) are practically equal to one for all 
planets. It agrees with Tycho Brahe astronomical observations used by Kepler. For a 
hundred years (415 "periods" of Mercury) the advance of Mercury's perihelion is nearly 
(1 - 71) • 360 • 415 • 3600" pa 7". 175. The relations (3.9), (3.22) imply the perihelion angle 

<t>k;i ~ fe; o + M 1 + 2 _1 ^a 2 c- 2 (l - e 2 )- 1 ), I = 0, ±1, ±2, .... (3.23) 

The substitution of the relations (3.14), (3.22) into the equality (3.9) yields 

e k cos ((l - 2- 1 o; 2 a 2 c- 2 (l - e 2 )" 1 ) (0 fc (t) - fc;O )) « a k (l - e 2 k )r^(t) - 1, (3.24) 

k = 1,...,9. In the general relativity the orbits ([4], Chap. 40, Sec. 40.5, relations (40.17), 
(40.18)) are the approximate solutions of the geodesic equation for the chosen metrics ([4], 
Chap. 40, Sec. 40.1, relation (40.3)). These orbits are the orbits (3.24) with the perihelion 
angles (f) k;0 = and with the precession coefficients 1 — 3w^a|c _2 (l — e|) _1 instead of the 
precession coefficients 1 — 2 -1 cj 2 .a 2 ,c~ 2 (l— e 2 ,) -1 . It seems that the perihelion angles are missed 
in [4]. We note that 3w 2 a 2 C - 2 (l - ej)- 1 ■ 360 • 415 • 3600" pa 6(1 - 71) • 360 • 415 • 3600" pa 
6 • 7". 175 = 43". 05. Does the orbit (3.24), k = 1, or the orbit ([4], Chap. 40, Sec. 40.5, 
relations (40.17), (40.18)) agree with the observed Mercury's orbit? From the paper ([8], 
p. 361) we know: "Observations of Mercury are among the most difficult in positional 
astronomy. They have to be made in the daytime, near noon, under unfavorable conditions 
of the atmosphere; and they are subject to large systematic and accidental errors arising 
both from this cause and from the shape of the visible disk of the planet. The planet's path 
in Newtonian space is not an ellipse but an exceedingly complicated space-curve due to the 
disturbing effects of all of the other planets. The calculation of this curve is a difficult and 
laborious task, and significantly different results have been obtained by different computers." 

Substituting the relations (3.20), (3.21) in the equality (3.13) and introducing the pa- 
rameter without physical measure we get 

Ofe V fe (r fc ) pa 1 + e fc sinTfc, 
u k t k (Tk) ~ n - r fe;0 - e fc (l - u 2 k a 2 k c~ 2 ) cosr fc , 

r fc = u k a k £ k (l + 2- 1 u 2 k a 2 k c' 2 ) £; = !,..., 9. (3.25) 
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If we neglect the values wf.af.c -2 , then the solutions (3.24), (3.25) of the equations (3.3) 
coincide with the solutions of the Kepler problem. Let us define the constant r k; o in the 
second equality (3.25) by choosing the initial time moment t k (0) = 0. Then the equalities 
(3.25) have the form 

a k lr 'k(n) ~ 1 + e fc sinr fc , 
0Jktk(n) w r k - e k (l - uj 2 k a 2 k c~ 2 ) (cos r k - 1) , fc = l,...,9. (3.26) 

Let the direction of the first axis be orthogonal to the vector M(xi). Let the direction of 
the third axis coincide with the direction of vector M(x 3 ). Then the second axis lies in the 
plane stretched on the vectors M(xi) and M(x 3 ). Due to the relations (3.6) 

x\{t) — ri(t) cos (j>i(t), x\(t) = — ri(t) cos Q\ sin <f>i(t), xf(t) — ri(t) sin 9\ sin (/>i(t), 

x\(t) = r 3 (t) cos0 3 (t), x 2 3 (t) = r 3 (t) sin0 3 (t), x 3 3 = (3.27) 

where the inclination of Mercury orbit plane 6\ = 7° and the values r k (t),<f) k (t), k = 1,3, 
satisfy the equations (3.24), (3.26). For the definition of Mercury and the Earth trajectories 
it is necessary to define the perihelion angles 0i ;O , 3; o in the equations (3.24). 

" Observations of Mercury do not give the absolute position of the planet in space but 
only the direction of a line from the planet to the observer." ([8], p. 363.) The advance of 
Mercury's perihelion is given by the angle 

(xi(ti(n,i)) -X3(t 3 (T3,i)),xi(ti(Ti i2 )) -x 3 (t 3 (r 3j2 ))) 



cos a 



-X3(t3(r 3) i))||xi(ti(r li2 )) -x 3 (i 3 (T3,2))r 
c(ts(T3,fc) -*i(ri )fc )) = |xi(ti(ri )fc )) -x 3 (t 3 (r 3ifc ))|, k — 1,2, 

ii(ri, 2 ) - *i(ri,i) < 100T 3 < t^) - ^(n.i) + T x (3.28) 

where the parameters t^i, ti j2 are defined by Mercury's perihelion points, the parameters 
73,1, r 3,2 are the solutions of the second equation (3.28), the numbers T±, T 3 are the orbit 
"periods" of Mercury and the Earth. The quotient T 3 /Ti of the Earth and Mercury orbit 
"periods" is approximately equal to 4.15. 

By making use of the equations (3.26) we obtain the parameters corresponding to Mer- 
cury's perihelion points: 

ar lr i( r i,fc) ~ 1 - ei, 
wifiKfc) w 7T (2l k + 3/2) + ei(l - cofalc' 2 ), 

n, fc w 7T (2l k + 3/2) , k = l,2, (3.29) 

where l k are the integers. The first relation (3.29) coincides with the equality (3.16). 

According to the book ([4], Chap. 25, Sec. 25.1, Appendix 25.1) c~V = 275.8-10~ 17 m- 1 , 
c^cos = 66.41 • 10- 17 m _1 , a x = 0.5791 • 10 n m, a 3 = 1.4960 • 10 n m. The substitution of the 
second equality (3.29) into the third relation (3.28) yields h — h — 415. 

Due to the second relation (3.28) 

X3(t 3 (T3,fc)) = x 3 (ti(r ljfc )) + c^lxi^i^fc)) - x 3 (t 3 (r 3 , fc ))|v 3 (t 3jfe ), k = 1,2. (3.30) 

The Earth speed is small compared with the speed of light: c _1 |v 3 | w c _1 w 3 a 3 w 0.9935-10 -4 . 
We neglect this value (arcsin 10" 4 w 0°.0057). Then the relations (3.28), (3.30) imply 

cog _ (xi(ti(n,i)) -x 3 (t 1 (r 1;1 )),x 1 (t 1 (r 1 , 2 )) -x 3 (t 1 (r 1 , 2 ))) 
~ |x 1 (t 1 (r lil ))-x 3 (t 1 (r 1 , 1 ))||x 1 (t 1 (r li2 ))-x 3 (t 1 (r li2 ))| 
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where the parameters Ti^, k — 1,2, are given by the third relation (3.29) and the relation 
h = h +415. 

Let us consider Mercury's perihelion points corresponding to the integers l\ = and 
l 2 = 415. The substitution of the values corresponding to the Mercury's perihelion, defined 
by the first equation (3.29), into the equation (3.24) yields 

afV (tt (2/ + 3/2)) w 1-ei, 
0i (*i (tt (21 + 3/2))) w 0i ;O + 2tt/ (l + 2~ 1 w 1 2 a 2 cr 2 (l - e 2 )" 1 ) (3.32) 

since the value io 2 a\c~ 2 m 2.5509 • 10~ 8 is negligible. We substitute the time, defined by the 
second relation (3.29), into the second relation (3.26) for the Earth 

a 3 lr 3(+3(0) ~ 1 + e 3 sin r 3 (Z), 
u 3 ti (tt (21 + 3/2)) « w 3 o;f 1 (tt (2/ + 3/2) + ei(l - u; 2 a 2 c -2 )) « 

r 3 (Z) - e 3 (l - w^c" 2 ) (cosr 3 (0 - 1) . (3.33) 



Solving the second equation (3.33) we get r 3 (0) « 1.1748, r 3 (415) « 629.09. Substituting 
these values in the first equation (3.33) we have a 3 _1 r 3 (r 3 (0)) ~ 1.0157, a 3 1 r 3 (r 3 (415)) ~ 
1.0118. We substitute the first equation (3.33) in the equation (3.24) for the Earth 

The function in the right - hand side of the equation (3.34) is monotonic with respect to 
the variable e 3 on the interval < e 3 < 1. Calculating the values of this function at the 
points e 3 = 0, 1 we get the estimation for the module of this function which implies that 
the equation (3.34) has a solution. Substituting the solutions t 3 (Z), I = 0,415, of the second 
equation (3.33) in the equation (3.34) we get the angles in radians 

3 (h (tt (2 • + 3/2))) w 3;O + 2.7521, 

3 (tt (tt (2 • 415 + 3/2))) w 3;O + 2.3544 + 2tt • 99 (l + 2- 1 w 2 a 2 (T 2 (l - e 2 )- 1 ) (3.35) 

since the value u 2 a\c~ 2 0.9870 ■ 10~ 8 is negligible. Substituting the radii and the angles 
(3.32), the radii (3.33) and the angles (3.35) in the equations (3.27), (3.31) we get the 
equation 

cos a(0, 415) w (((1 - ei) cos 0i ;O - 1.0157a 3 a^ cos(0 3;O + 2.7521)) 

x ((1 - ei) cos(0i ;O + 4157rw 2 a 2 £T 2 (l - e 2 ) -1 ) 
-l.OHSaga^ 1 cos(0 3;O + 2.3544 + 997rw 2 a 2 cT 2 (l - e 2 )" 1 )) 
+ (0.99255(1 — ei) sin 0i ;O + 1.0157a 3 a7 1 sin(0 3;O + 2.7521)) 
x (0.99255(1 - ei) sin(0 1;O + 415vrw 2 a 2 c" 2 (l - e 2 )" 1 ) 
+(1.0118a 3 a7 1 sin(0 3;O + 2.3544 + 99vrw 2 a 2 cr 2 (l - e 2 )" 1 )) 
+0.01485(1 - e x ) 2 sin 1;O sin(0 1;O + 4157rw 2 a 2 c- 2 (l - ej)' 1 ))) 
x (((1 - ei) cos0i ;O - 1.0157a 3 a7 1 cos(0 3;O + 2.7521)) 2 
+ (0.99255(1 - ex) sin0 1;O + 1.0157a 3 ar 1 sin(0 3;O + 2.7521)) 2 
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+0.01485(1 - d) 2 sin 2 0i ; o) _1/2 
x (((1 - ei) cos(0i ;O + 4157rw 2 a 2 (T 2 (l - ejy 1 ) 
-l.OHSaaar 1 cos(0 3; o + 2.3544 + 997rw 2 a 2 c" 2 (l - e 2 )" 1 )) 2 
+ (0.99255(1 - e ± ) sin(0 1;O + 4157rw 2 a 2 c~ 2 (l - ejy 1 ) 
+ 1.0118a 3 a^ 1 sin(0 3;O + 2.3544 + 99tt + w 2 a 2 (T 2 (l - e 2 )" 1 )) 2 
+0.01485(1 - d) 2 sin 2 (0i ;O + 4157ra; 2 a 2 c" 2 (l - ejy 1 ))- 1 / 2 . (3.36) 

The perihelion angles 1;O , 03 ; o are needed. For the orbits ([4], Chap. 40, Sec. 40.5, relations 
(40.17), (40.18)) the perihelion angles (fik-o = 0. Let the perihelion angles 0i ;O , 03 ; o in the 
equation (3.36) be equal to zero. Then a(0, 415) = 17°. 889. According to the book ([4], 
Chap. 40, Sec. 40.5, Appendix 40.3), the observed advance of Mercury's perihelion is 
1°. 55548 ± 0°. 00011 for a hundred years. The angle a(0,415) = 17°. 889 is not small. In our 
opinion for the experimental verification of the general relativity it is necessary to obtain the 
advance of Mercury's perihelion, observed from the Earth, by making use of the Mercury 
and Earth orbits ([4], Chap. 40, Sec. 40.5, relations (40.17), (40.18)) calculated without 
Newton gravity theory. The orbits ([4], Chap. 40, Sec. 40.5, relations (40.17), (40.18)) 
are the orbits (3.24) with the perihelion angles <pk;o = and with the precession coefficients 
1 — 3w 2 a 2 c~ 2 (l— e 2 )" 1 instead of the precession coefficients 1 — 2~ 1 w 2 a 2 ,c~ 2 (l — e 2 ,) -1 . It seems 
that the perihelion angles are missed in [4]. In order to calculate the advance of Mercury's 
perihelion we need to know also the time dependence of the orbit ([4], Chap. 40, Sec. 40.5, 
relations (40.17), (40.18)) radius. 
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